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Y. C. Zou, Z. G. Dai and D. Xu 1 



ABSTRACT 

By considering synchrotron radiative process in the internal shock model and 
assuming that all internal shocks are nearly equally energetic, we analyze the 
gamma-ray burst (GRB) emission at different radii corresponding to different 
observed times. We apply this model to GRB 050904 and find that our analytical 
results can provide a natural explanation for the multi-band observations of GRB 
050904. This suggests that the X-ray flare emission and the optical emission 
of this burst could have originated from internal shocks being due to collisions 
among nearly-equally-energetic shells ejected from the central engine. Thus GRB 
050904 appears to be a burst with super-long central engine activity. 

Subject headings: gamma rays: bursts - hydrodynamics - relativity - shock waves 



Introduction 



The gamma-ray burst (GRB) 050904 was an explos ive event at redshif t of z = 6.29±0.01 



which has been m easured through various methods (IKawai et al. 



Haislip et al.ll2005l ). After the Swift trigger (jCummings et al 



2005; Price et al. 2005 



20051 ). multi- wavelength ob- 



servations of this high-redshift burst were performed. The 7-ray iso tropic-equivalent energ y 



of this bur st E^i„ n was between 6.6 x 10 53 ergs and 3.2 x 10 54 ergs (jCusumano et al.ll2005l ). 



Boer et al.l (120051 ) observed the optical emission at times of 86 s to 1666 s. The optical a f- 
terglow and its spectrum were analyzed quickly (IHaislip et al.ll2005l ; iTagliaferri et al.ll2005l ). 
A break of the optica l afterglow light curve was also observed at about 2.6 ±1.0 days 
(ITagliaferri et al.ll2005l ). The optical emission appears to be understood in ei ther both the 
late internal shock model or the reverse-forward shock model JWei et al.lhoOon . 



One of the most remarkable features of this b urst is a long l a sting varia bility of the X-ray 
emission, showing several X-ray flares (XRFs). I Watson et al.l (120051 ) and ICusumano et al. 
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( 120051 ) suggested that this variability may be due to a l ong-lasting activity of the central 



engin e. The XRFs we r e also found in some othe r GRBs (INousek et al.l 12006c lO'Brien et al. 



2006|). Burrows et al.l (|2005h . Fan fc Weil (12005( 1. IZhang et all (120051 ) and 



Wu et al 



" fEooih 

assumed late central engine activities to interpret the XRFs. The plausible origin models 
of XRFs were re cently proposed, e.g., the magnetic activity of a newborn millisecond pulsar 
( Dai et al.ll200q). f ragmentation of a neutron star by a bla ck hole in a compa ct object binary 
( Faber et al. 2006 ). fragmentation of an accret ion disk (Perna et alibood ). and magnetic 



barrier-driven modulation of an accretion disk (jProga et al.ll2006l ). However, there has not 



yet been any explicit and detailed investigation on this highly-varia ble light curve. Enlight- 
ened by the internal shock model to interpret the 7 -ray emission (jRees fc Meszaros 1994 ) 
and prompt optical emission (IMeszaros &: Reedll999l ). we assume that all the highly- variable 
X-ray emission (even at very late times) originate from internal shocks. These shocks oc- 
cur when many faster shells with equal energy and equal mass catch up with a slower shell 
with much more kinetic energy and mass. By fitting the peaks of the X-ray and optical 
light curves, we find that these internal shocks can lead to the observed X-ray and optical 
emission. We describe an internal-shock model in $2j fit the peaks of X-ray and optical light 
curves of GRB 050904 in SJ31 We summarize our results in £ffl 



2. Model 



We assume two cold shells labeled with 1 and 4, whose isotropic kinetic energies are 
Ei and E±, widths are A^o and A^o in the cosmological rest frame, and Lorentz factors 
are 71 and 74 (74 3> 71 3> 1), respectively. The shells come from the central engine with 
time lag T (also in the rest frame). Shell 4 catches up with shell 1 at radius r int ~ ^fcT 
and then two internal shocks occur, that is, a forward shock that propagates into shell 1 
and a reverse shock that propagates into shell 4. Four regions are divided by the reverse 
and forward shocks: unshocked shell 1 (region 1), shocked gas of shell 1 (region 2), shocked 
gas of shell 4 (region 3) and unshocked shell 4 (region 4), with a contact discontinuity 
(CD) surface separating region 2 and region 3. Figure [TJ gives a sketch of this model. The 
number densities of region 1 and region 4 at radius r are n\ = i?i/(47rr 2 7^Ai i om p c 2 ) and n± = 
E±j (4:7ir 2 jl A 4i om p c 2 ) respectively. To produce plenty of X-rays and 7-rays, it is required that 
the reverse-shocked materials are relativistic, i.e., 7 34 ~ -^(-^Y^ 2 /~ 1//4 3> 1 ( Sari fc Piran 
19951 : IWu. Dai fc Huangjl2003l : IZou. Wu fc Daill2005l ). where 734 is the relative Lorentz factor 
between regions 4 and 3, and the number density ratio is defined as / = 714/ni. 

T o calculate the w i dth of the coasting shell, we define the spreading radius -Ra = 
7 2 A ( IMeszaros fc Reed Il993l ; iPiran. Shemi &: Narayanl Il993l ). For the internal shocks of 
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our interests, if r < min(i?A,i = 1.0 x 10 15 7 2 1 5 Ai i o,i2cm, Ra,a = 1-0 x 10 16 7f gA^o.iocm), 
the width of the shells can be considered to be constant, and for r > max(i£A,i, Ra,a)j the 
width should be regarded as r/7 2 and r/7 2 for the spreading of the shells. The conventional 
notation Q k = Q/10 k with cgs units is used in this paper except for special explanations. 
We now consider these two cases. 



2.1. No-spreading case 

At radius r < min(i?A,i, Ra,a), for a certain set of parameters Ai )0 ,i 2 = A 4i o,io = #1,54 = 
■#4,51 = 71,1.5 = 74,3 = 1, one obtains 7^/ ~ 2.5 x 10 -2 ^ 5 1 4 E 4)51 A 1i q )12 A^J j10 < 1 (where 
714 is the relative Lorentz factor between regions 4 and 1), implying that the forward shock 
is Newtonian, and 7i 4 // — 2.5 x 10 6 E^, M E^^ A7n i9A 4 .n.in 7i~i 4 ^7^ ^> 1, which shows that 



the reverse shock is relativistic. Using the shock conditions ( jBlandford fc McKedll976l ). and 



keeping equali ty of pressures and Lorentz factors of the materials at two sides of the contact 



discontinuity (ISari fc Piranl Il9951 ). we can obtain the number densities n 2 ,n 3 , and energy 



densities e 2 , e 3 of the shocked regions in the comoving frame, 

n 2 ~ 1.0xl0 10 (l + 2 ) 2 E li5 4A7 i i 12 7!: 1 6 .5^ 2 2 cm- 3 , (1) 
n 3 ~ 9.3 x 10 6 (1 + ^^^siA^o 7^.5 7 4 r 3 1 i- 2 2 cm" 3 , (2) 
e 2 = e3~2.2xl0 5 (l + z) 2 E 4i 5iA 4i0il0 7j 5 t- 2 2 ergcm- 3 . (3) 

The Lorentz factor 7 2 = 73 ~ 71 ~ 31.671,1.5. Here is the observer's time. 

Assuming that the reverse shock lasts from r int to r int + Sr, we find 5r ~ 1.0 x 
10 13 7i 2 i. 5 A 4 , 10 cm, which corresponds to a pulse from the beginning to peak. In the ob- 
server's frame (dt® = (1 + z)dr/(2 , y 2 c)), this lasting time is 

5t® ~ 0.17(1 + 0)A 4A ioS. (4) 

By considering the synchrotron emission at the reverse shock-crossing time, which cor- 
responds to the peak emission time, we obtain the characteristic frequencies ^ m ,2, P and 
v m ,3, P , cooling frequencies v c ,2, P and u c ,^,pi synchrotron self-absorption frequencies v a ,2,p anc ^ 
Vn^p, and the maxim um flux densities f v ,max,2, P and f u ,max,3, P in the observer's frame (see 



Zou. Wu &: Daill2005l . for original formulae of synchrotron emission) 



5 51 

j 1 2 1 — a 2 a — 2 2 2 1. / \ 

v m,2,p — 1.3x10 #1,54 #4 2 ,51 ^1,0,12 ^4,0,10 e B,-l e e, -0.5 7l, 1.5 Cl/6 ^e, 2 Hz, (5) 



v m , 3 ,p ^ 1.5 x 10 17 E\ M A 4 | 10 e\_ x e 2 _ . 5 7 J. 5 74 2 3 tf /6 Hz, (6) 
v c> 2 >P = v^p - 1.2 x 10 15 (1 + z)- 4 A~J 10 e-Li 7i, 1.5 4,2 Hz, (7) 
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v aXv ~ 3.2 x 10 i4 (1 + z)-ts EJt £| 51 A- il2 A 4) ° 10 ej.^ 7l /. 5 C* /6 f~ ? Hz{8) 



v, 



a,3,p 



,14 1 - ' 



1 



1.2 x 10 14 (l + ^)-3£| 51 A^| )10 7x7.5 7 4 7 VI Hz, (9) 

/„,ma X ,2, P ^ 8.0 X 10 2 L> 2 - 8 2 4,54 #4,51 (1 + Z? A"^ 7^.5 ^2 ty, (10) 

/,,ma X ,3, P ^ 0.47 (1 + ^) 2 I>m 2 451 A W0 4,-1 Tu.6 7^3 ^2 J Y> 

where Ci/6 = 6(p — 2)/(p — 1), p is the spectral index of the shock-accelerated electrons, -D 
is the luminosity distance from the burst to the observer (as a function of redshift z), and 
e e and €b are the fractions of the internal energy density that were carried by electrons and 
magnetic fields, respectively. In this paper, p is taken to be 2.2, and e e and €b are equal for 
both region 2 and region 3. Note that the expressions of ^ a ,2,p and v a ,3,p are only valid in the 
cases of v Cj2 , p < v m ,2, P < v a,2, P and v aAp < u cAp < v m ^ v . 

From the above equations, we see v c ,2, P < v m ,2, P < v a,2,p and v a ,3, P < v c,3, P < ^m,3,p for 
typical parameters. The flux density in the fast-cooling case at the peak time is then given 
by 

2.4 x 10 s (1 + z)- 1 D 28 2 E^ €^_ x 7 8 )L5 v\ 8 t% 2 Jy, v < v c , 2 , P < v m ,2, P < v a ,2, P , 

7.1 XlO 9 (1 + ^28^5^0,^ ^ Ve,2, P < V < Va,2, P , (12) 

2.0 x 10- 3 D 2S 2 E~i Ej M A^ 12 A"| 10 ef_, ej_ 0>6 7^.5 "ifi C 1/6 %f Jy, ^,2, P < ^ m ,2, P < ^a,2, P < " 
for the forward shock emission, and 



3,p 



) 1 D 2 i E A 51 e B L _ 1 t| 1>5 z/ 2 8 4 j2 Jy, 
10 2 — ~ ~ ~~ ~ 4 - - 

) 3 -^28 ^4,51 A 4,0,10 e B,-l 7l,1.5 ^18 ^©,2 JY: 

3 _3 _1_ _1 1 

TPA A 4 - 4 -.2 ..— 1 .. 2 J-2 T,. 



3.3 x 10 8 (1 + 2) 1 D 2 £ EA X e B x _ x 7 8 ;1 5 v\ 8 4,2 Jy, v < "c,3, P < v a ,3, P < v m ,3, P , 

i 1' i k r> 



9.6 x 10 9 (1 + z)t D 28 2 E A ^ A| 010 7^.5 1/& t^ 2 Jy, z/ Cj3iP < v < u a ^ p < u mAp , 

(16) 

< V < V m> 3, P , 



3 _ 3 _i_ ' I ' 

1.6 X 10~ 2 DnTg 2 E 4 3 51 A 4 £ 10 e fl 4 _ 1 7?, 1.5 74> *V 4,2 JY, ^,3,p < ^a,3,p 

2 21 _21 J_ 6 _2 1 6 _11 _J_ 

5.2 x 10 _ £> 28 ^4,51 A 4i o > ° 10 e^ _ 1 e| j _ 5 7 1) i 5 7 4 5 j 3 Ci/6^18 10 *©,2° Jy, ^,3^ < ^a,3, P < ^m,3, P 
for the reverse shock emission. 



2.2. Spreading case 

In the case of r > max(i?A,i, -Ra,4), the widths of the shells become A x ~ r/7 2 and 
A 4 ~ r/74 because of the spreading effect. For the same parameters as §2.1[ 7 2 4 / ~ 



< v. 
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O.S-E-l 54-E , 4,5i7 1 1 5 74 3 < 1, showing that the forward shock can be considered as Newto- 
nian approximately, and 7i 4 // — 2.5 x lO^i^i?^ 77/1.574 3 3> 1, implying that the reverse 



shock is still relativistic. 

Using the shock conditions as in §2.14 we obtain 

n 2 ~ 1.7 x 10 9 (l + z) 3 E li54 7i;i.5^ 3 2cm- 3 , (14) 

n 3 ~ 1.5 x 10 7 (1 + zf £4,51 7J.5 74,3 ^e 3 2 cm " 3 > ( 15 ) 

e 2 = e 3 ~3.7x 10 5 (1 + E 4 , 5 i 7J.5 7 4 2 3 ^2 erg cm" 3 , (16) 

72 = 73 — 7i — 31.671,1.5. (17) 



The distance for the reverse shock to cross shell 4 is Sr ~ 6.0 x 10 12 (1 + z) X 7 4 x 5 t ej2 cm, 



and the rising time of the pulse in the observer's frame is 8t® ~ 0.27 43 t ej2 s. 
The corresponding emission values at the peak time are given by 

v mXv ~ 1.6xl0 13 (l + ^)^rj 4 4,5i4,-i<-o.5 7j.5 74%Ci%^! Hz ; (18) 

v m>3 , p ~ 2.0xl0 17 (l + z)^| 51 e| i _ 1 e e %. 57l 7. 5 74 3 3Ci%vl Hz , (19) 

3 3 5 

^ C ,2, P = Vc^ v ^ 1.5 x 10 15 (1 + z)-i E^l e~ B *_ x 7l 7 iL5 74 , 3 1| >2 Hz, (20) 

v aXv ~ 6.4 xl0 14 (l + z)- ^,5i^i<-o.5 7r:i 2 .5 74:3Q 6 V2 Hz, (21) 

^, 3lP ~ 2.0xl0 13 (l + 2)^45i4,-i7i:i.5 74"l^jHz, (22) 

/,,max,2, P ^ 0.033(l + z)lD 2 - 8 2 ^ 54 E 4i5 i e | i „ l7 r >1 2 5 te!jy, (23) 

/,, max ,3, P ^ 0.6 (1 + ^)t^ 2 - 8 2 4,51 4,-i 7i:i. 5 v!jy- (24) 



Note that the expressions of ^ a ,2, P and is a ,3, P are valid only in the cases of z/ mj2iP < is c ,2, P < v a,2, P 
and v a ^ p < j/ C)3)P < f mj3jP respectively. 

In the case of v m ,2,p < v c,i,p < ^0,2^5 the corresponding flux density of the forward shock 
emission is 

{2.1 X 10 7 (1 + Z) ^54 £4,51 e e -0.5 71,1.5 74,3 ^18 Cl/6 t®,2 J y> ^ < ^m^p < f C ,2,jp < ^a,2,p, 

5.2 X 10 9 (1 + z) f D 2 ~ 8 2 e~ ^ 7i f ! 5 7 4 ~f 4 *| 2 Jy> "mAp < ^ < ^a,2, P , (25) 

21 -JL 1 _L £ _22 2 _ii £ _23 

1.7 x 10- 2 (1 + 2)55 D 28 2 E x $ E* S1 eg^ e*_ . B 7i jl 10 5 7|, 3 "is 10 C 1/6 *e,!° Jy, "m^p < v c ,2, P < »a,2, P < 

As the cooling frequency v c ,2, P exceeds the self- absorption frequency v a ,2, P at time £® ~ 70(1 + 

\ 7-1—7/134 t-i37/67 a 7/134 a —1/67 57/134 6/67 -160/67 A 6/67 j.' f , 

^^1,54 ^4,5i A i,o,i2 A 4,o,io e J 3,-i e e,-o.57i,i. 5 Ci/ 6 S ' the corresponding emission after this 
time is then 

fv,2,P - 



(-e -0.5 71,1-5 74,3 ^18 Cl/6 *©,2 ^ V < ^m^.p < ^,2^ < ^,2^, 

3 _I 2 _I 5 II 

7.1 x 10 9 (1 + z)f D 28 2 J B 4 5 4 1 7^ 5 7 4 3 2 z/^ t * a Jy, z; m ,2, P < v < v aXv < u Cj2 , P , 

14 r,_X'5 i 6 _31 _36_12 

4.4 x KT 1 (1 + z)t D 28 E^l El bl e% x e e 5 _ 5 7lil 5 5 7 | )3 V Q/& *©,2 J y> "m,2,p < ^a,2, P < ^ < f c ,2, P , 
21 —X 1 _L 5 _22 2 _11 6 _23 

[ 1.7 x 10- 2 (l + 2)*> D 28 2 £ 4 51 e |»_ 1 e J_ 05 7lj w 7 | 3 ^10 C K /6 Jy> , m ^ p < ^ < ^ < „. 

The flux density of the reverse shock emission in the case of v a ,z,p < u c,3, P < v m,3,p is 
given by 

fv,3,p — 

3.5 x 10 s (1 + z)- 1 D 2 2 E^h e B 1 _ 1 jf >15 uf 8 i|, 2 Jy, v < v a ,3, P < ^ C ,3, P < v m ,3, P , 

11 —Ii _I I _2 

5.3 (l+z)s- D 28 2 Ej 51 e B ,-i 7ii 3 5 7 4 ,3 3 U 18 *©,2 Jy> ^,3, P < f < ^c,3, P < fm,3,p, . , 

3 _ 2 s -i 1 I _i _i (27) 
2.3 x 1CT 2 (l + 2)i D 2S Ef hl e B 4 _i 7i,i. 5 7 4 2 ,3 ^is' *©,2 J y> ^,3, P < ^c,3, P < ^ < fm,3, P , 

21 21_L 2 _5236 _11_23 

8.9 x 10- 3 (1 + 2)20 D 2 2 E% x e%_ x ej_ 0>5 7J.5 7$ C x % V° *©,f J y> ^3* < ^,3, P < ^ m ,3, P < v. 

At time £ ffi ~ 3.4 x 10 2 (1 + ^E^ 2 eB^ 1 eef-o.57u. 9 5 /8 74,3 2 Ci 1 /6 s > "c,3, P = ^m,3, P - After that time, 
the order of the break frequencies becomes ^ a ,3, P < v m,3,p < v c,3, P - In this case, the flux 
density becomes 

fv,3,p — 



3.4 x 10 9 (1 + z) D 28 e e _ .5 71,1.5 74,3 Ci/6 *T 8 *©,2 J y> ^ < v a,3j> < U c,3,p < Vm,3,p, 

4 1 2 4 2 1 



9 (1 + Z) D 28 e e> _ .5 7? 1.5 74,3 Cl/6 ^1 2 8 % 

7 4 1 _2 __4 _2 j_ 

1.0 (1 + z) 3 D^ 2 £f 51 ^ e e I 5 7l)1 5 5 7 7j C v i ^i 3 8 ^,2 J y> ^, 3 , P < v < v cAv < z/ m , 3 

14 „ 2 i 6 „ 2 2 _3 _12 

2.3 x KT 1 (1 + z)t £>" E 4 5 51 e| x e|_ 5 7^ 5 7^ Cj 6 V t e j Jy, ^.,3, P < ^ C ,3, P < v < v m ,3, P , 

21 21 J_ 6 _ 5 23 _11 6 _23 

8.9 x 10- 3 (1 + z) 20 D" 2 eg^ ej_ . 5 7i,L 74 1 ,° 3 ^i 8 10 C x % ^©,2 Jy> ^, 3 , P < ^ C ,3, P < ^ m ,3, P < 

From equations (1271) and ( 1281) . we can see that the flux densities are the same for the case 
in which the observed frequency exceeds the three break frequencies. 

We consider that many shells like the foregoing shell 4 collide with shell 1 at different 
radii r int , where r int indicates the observed time by the relation £ e ~ (1 + z)r/(2 7 |c). Thus 
the peaks of emission obey a power law temporal profile. Under more realistic conditions, 
"shell 1" may be sole, proceeding at the front, and many shells like "shell 4" ejected from 
the central engine continually overtake "shell 1" (see Figure UJ). As the energy of "shell 1" 
is assumed to be much greater than that of "shell 4", the dynamics of "shell 1" after the 
collision can be almost unchanged, and thus the above relations are still applicable. In the 
popular collapsar model, "shell 1" can be understood as due to the envelope of a collapsing 
massive star. Because a relativistic outflow needs to accumulate enough energy to break out 
of the envelope and is inevitably polluted by the baryons in the envelope, this outflow (i.e., 
"shell 1" ) should be more energetic but have less Lorentz factor comparing with "shell 4" . 
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GRB 050904 



Using the above model, we fit the peaks of the X-ray light curves of GRB 050904 in the 
observer's frame. We consider the model parameters E\ = 10 55 erg, E± = 1.1 x 10 51 erg, 71 = 
31.6, 74 = 1455, A 1>Q = 10 12 cm, A 4i0 = 2.0 x 10 8 cm, e B = 0.1, e e = 0.316, and p = 2.2 and the 
cosmology with Vt m = 0.23, £l\ = 0.73 and H = 71 kms -1 Mpc -1 . These parameters lead 
to a Newtonian forward shock and a relativistic reverse shock. As shown in Figure [21 the 
solid line with four segments is fitted analytically. The flux is an integral of flux density in a 
frequency range of 4 x 10 16 Hz to 2 x 10 18 Hz, which corresponds to XRT's 0.2-10 keV band. 
The temporal breaks are due to crossing of the three break frequencies. 

The first segment represents the emission from a relativistic reverse shock with a con- 
stant shell width, which satisfies equation (|T3|) with temporal index 1/2. At time ~ 42 s, the 
corresponding radius is R(42 s) ~ 2j^ct & /(l+z) ~ 3.5 x 10 14 cm, while the spreading radii of 
region 1 and region 4 are R s i ~ 7i A 10 — 10 15 cm, and R s ^ ~ 4.2 x 10 14 cm respectively. The 
three radii are approximately equal, and thus R(A2 s) can be taken as the common spreading 
radius of regions 1 and 4. 

After the first break time, regions 1 and 4 should be spreading and the flux density 
is described by the third expression of equation (T2T|) . The temporal index is —1/4. Up to 
~ 350 s, the frequency u m ^ p ~ 2.6 x 10 17 Hz decreases to the observed frequencies of XRT, 
and then the order of the break frequencies becomes f ,3,p < v c,3, P < ^m,3,p < v or v a ^ yP < 
v m ,3 lP < v c,3, P < v. The flux density is given by the last expression of equation fl27|) . and the 
temporal index of flux is —23/20. After ~3x 10 4 s, the order is v a ,3, P < v m,3, P <v< ^ C ,3, P , 
and the temporal index becomes —12/5. 

Observationally it was reported that the photon indices are ~ — 1.3 and ~ — 1.9 (whose 
correspond ing spectral indices are —0.3 and —0.9 resp ectively) at early and late times re- 



spectively (ICummings et al.l 120051 ; IWatson et al.l 120051 ). The observer's time at which the 



spectral index changed is about 350 s. This time can be naturally understood as the crossing 
time of v m to v. Before this crossing time, v c < v < v m and the spectral index is —1/2, and 
after this crossing time, v c < v m < v and the index is —1.1. Figure [3] shows the spectral 
index evolution and a sketch for the evolution of spectral energy distribution. These spectral 
indices are generally slightly less than the observed ones, because our model only gives the 
spectral indices at the peak times but the observations were performed in relatively broad 
ranges of time. At the be ginning of each pair of reverse and forward shocks, the cooling 



frequency v c is very large (jZou. Wu fc Daill2005l ). because this frequency is proportional to 



t~ Q 2 , where t co is the comoving time of the shocks. Thus v c can be larger than v during the 
initial stage of the shocks. In the case of v < (u m , u c ), the spectral index is 1/3, which makes 
the observed spectral index larger than the expected values —0.5 and —1.1. 
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The optical data during the same period were reported (jHaislip et al.ll2005t iBoer et al. 



20051 ). The observed J-band (2.47 x 10 14 Hz) data are plotted in Figure [2] in units of magni- 
tude. Using the above parameters, we plot a line representing the J-band peak flux density. 
However, the J-band emission may be mainly contributed by the forward shocks. Since the 
forward shocks are Newtonian, the emission contributes more at optical band than does at 
the X-ray band. Furthermore, as the number density of the forward-shocked region is much 
larger than that of the reverse shocks (see equations (IT4"|) and (TT5T) ) . the contribution to the 
J-band emission exceeds that of the reverse shocks. To show the contribution of the for- 
ward shocks, Figure H] duplicates Figure [2] but counts in the forward shocks. It is reasonable 
to see that the fitted peaks are greater than the observed data. As the fitting lines stand 
for the peak fluxes of the emission of internal shocks, these fluxes should be greater than 
the time-averaged observed data especially for the optical emission, which lasts a relatively 
longer time for the observation. If there are not reverse-forward shocks at a certain time, 
the observed flux density should be much less than the expected peak one. Therefore, the 
fitting solid line is basically consistent with the observed data. 

For the highest (fourth) J-band point (449s-589s) and the X-ray peak nearly at the 
same time, we suggest that they could arise from the same internal shocks, and the X-ray 
peak could occur only at time slightly earlier than ~ 440 s (we circled them in Figure H] for 
emphasis). There is a time lag (~ — 140s) betwee n the X-ray emission and the J- band 
emission. This may be caused by a spectral time lag (INorris. Marani fc Bonnellll2000l ). the 
different refractivities of photons in different energy bands during the propagation from the 
source to the earth, or a quantu m gravity effect (which is also d ue to different refractivities 
in vacuum for different photons, lAmelino-Camelia et al.l (119981) ). There ha s also been an 
example with a long time lag, GRB 980425 (INorris. Marani fc Bonnellll2000l ). which shows 
the probability of time lag with tens of seconds. 

We further show that our model is consistent with the observed data at some other as- 
pects. Firstly, the Lorentz factor of region 1 should not change significantly during the period 
of X-ray emission. Because faster ejected shells (region 4) have much less kinetic energy than 
the front slower one (region 1) and the forward shocks are Newtonian, the Lorentz factor and 
the particle number of region 1 do not increase significantly when it is overtaken by faster 
shells. On the other hand, we estimate the deceleration radius r d = [3E / ' {A-K^nrripC 2 )] 1 ^ ', 
where Eq is the isotropic-equivalent kinetic energy, and n is the medium density. Letting 
the gamma-ray isotrop ic equivalent energy be E lt i SO = 2.5 x 10 54 erg, which is consistent with 
the observation range (ICusumano et al.l 120051 ). and assuming the radiative efficiency r\ = 0.2 
and the medium density n = 3.16cm~ 3 , we obtain E = E 1)iso (l — rj)/r] = 1 x 10 55 erg. 
This energy is consistent with the parameters assumed above. And then the decelera- 

! cm, which corresponds to the observed time 



17pV3 -2/3-1/3 



tion radius is ~ 8.0 x 10 7 E 55 j 1 



5 n 0.5 
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~ 9.7 x 10 4 [(1 +z)/7.29]£/o557i.5^ 3n o.5 s - This deceleration time is larger than the X-ray 
duration time, indicating that region 1 is not decelerated by the ambient medium during 
this period. 



Secondly, the Lorentz factor at the jet-break time (ITagliaferri et al.l 120051 ) is about 
1/6 j ~ 17.5 ( 8j ~ 0.057). It should be greater at an earlier time. This is in agreement 
with our used parameter (71 = 31.6). 

Thirdly, we note the photon indices at the second break of the fitting X-ray peak line 
in Figure EJ The indices are —1.22 ± 0.10 during the time period of (149s, 223s), —1.5 ± 0.3 
during the time period of (223s, 303s) in the BAT data, and —1.13 ± 0.07 during the time 
period of (169s, 209s), -1.31 ± 0.06 during the time period of (20 9s, 269s), -1.34 ± 0.06 



during the time period of (269s, 371s) in the XRT data (Table 1 of ICusumano et al.l 12005 
where the time is divided by (1+z) to give the local observer's frame). It can be seen that the 
photon indices of the BAT data are greater than those of the XRT data. In our model, the 
frequency v m is crossing the observed frequency during this period. For the frequency orders 
u c < v < u m and u c < u m < u, the flux densities are scaled as f u oc z/~ 1//2 and f u oc z/~ n / 10 
respectively. Therefore, the spectrum is harder for the lower-frequency emission than for the 
higher-frequency emission, which is consistent with the observed data. 

Fourthly, the mean intrinsic column density is Nn = (2.30 ± 0.50) x 10 22 cm -2 , as 



measured by XRT (ICusumano et al.l 120051 ) . This high density cannot be explained by the 
ambient medium with assumed number density 3cm -3 . Because the radiation originates 
from the reverse-forward shocks and should pass through "shell 1" , we find that the column 
density of "shell 1" is Nn = E\/ (7i47rr 2 m p c 2 ), which is 1.67 x 10 23 cm -2 at radius r = 10 16 cm 
and 1.67 x 10 21 cm -2 at radius r = 10 17 cm. These values are basically consistent with the 
measured mean column density. Therefore, the materials in "shell 1" provide an explanation 
for the high column density measured by XRT. 

Finally, there are several low peaks in figure [2] at about 10 3 s. They may be caused by 
non-uniform Lorentz factors or energies of the ejected rapid shells. 

In short, the X-ray flares are understood as due to the internal shocks. These shocks are 
caused by the continually ejected faster shells that overtake the sole slower one with much 
larger kinetic energy and much more baryons. 
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4. Conclusions and discussion 

In this paper, we have proposed a simple internal shock model: a slower massive shell 
propagates in the front, and many faster shells with similar mass and energy catch up with 
the front slower shell at different radii. These collisions lead to reverse-forward shocks. The 
shocks produce the 7-ray and lower-energy emission by synchrotron radiation. We obtained 
a "normal" set of parameters to fit the peak fluxes of XRFs of GRB 050904, and found that 
our model is well consistent with the observations including the early optical data. 

For GRB 050904, we also found that the contribution of the forward shock emission 
is insignificant in the X-ray band, but significant in the J-band at early times, because the 
frequency of emission from the Newtonian forward shocks mainly lies in the optical band. 
However, we should note that, for other bursts, forward shocks may also be important for 
the X-ray emission. If the leading shell ("shell 1") is not more energetic than the followed 
faster shells ("shells 4"), the velocity of "shell 1" cannot be regarded as unchanged after 
some collisions and the simple model in section [2] is invalid. If the parameters differ, the 
deceleration radius may be smaller, and the external shock may become important. All the 
effects would change the profile of the XRF peaks. 

In general, we argued that some long bursts should have the same properties: the peaks 
are nearly constant before the time at which v m is just equal to the observed frequency u, 
and decline rapidly with temporal index about — (3p — 2)/4 after this time. Therefore, the 
7-ray emission of these bursts persists up to the time when u m = u, and disappears because 
of the rapid decline. In this case, the 7-ray emission ceases when the peak energy E p ~ hv. 
There are two other reasons for the disappearance of 7-rays. One is that the flux of 7-ray 
emission decreases below the detectable limits of telescopes when the frequency order is still 
v c < v < v m . The other is the turnoff of the central engine when v m is still larger than v. 
Only in the latter case, it can be said that the 7-rays duration is the burst duration. 

We would like to thank the referees for valuable comments that have allowed us to 
improve our manuscript and X. F. Wu, L. Shao, Y. F. Huang, B. Zhang and Y. Z. Fan 
for helpful discussions. The formulae were generated by the help of Maxima. This work 
was supported by the National Natural Science Foundation of China (grants 10233010 and 
10221001). 
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contact 
discontinuity 




shell 1 



Fig. 1. — Sketch of our model: A slower shell labeled "1" with much more energy expands 
in the front. Many similar faster shells labeled "4" with less energy and less number density 
catch up with "shell 1", and reverse-forward shocks occur. The gamma rays and early X-rays 
are mainly emitted from the reverse shocks. 
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Fig. 2. — Fitting to the peaks of X-ray and optical light curves of GRB 050904 (soli d line) . 
The observed X-ray data (points with error bars) are taken from ICusumano et al.l (120051 ) 
where the flux is divided by (1 + z) 2 and the ti me is multip l ied b y (1 + z) to give the 
observer's frame. The optical data are taken from iBoer et al.l J2005h (triangles with error 
bars, extrapolated to J-band by t he relation F„ oc v 1 ^ , which is the synchrotron emission 
in the case v a < v < v c < u m ) and lHaislip et al. J2005h (squares). The top panel shows the 
break frequencies as functions of time. Vertical dashed lines represent the crossing times of 
any two frequencies. 
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Fig. 3. — Spectral index evolution. The thick solid line is the expected spectral index at the 
peak of synchrotron radiation. The corresponding spectral energy distribution for the two 
index values (—0.5 and —1.1) are plotted schematical ly in the two sub-figures. The points 
with error bars are taken from ICusumano et al.l (120051 ) . 
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Fig. 4. — Same as Figure [2] except for counting in the contribution of the forward shocks. 
The solid line indicates the total peak flux density and the others are annotated in this 
figure. 



